Incommensurate spin Luttinger liquid phase in a model for the spin-Peierls materials 

TiOBr and TiOCl 
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In the present work we aim to characterize the lattice configurations and the magnetic behavior in 
the incommensurate phase of spin-Peierls systems. This phase emerges when the magnetic exchange 
interaction is coupled to the distortions of an underlying triangular lattice and has its experimental 
realization in the quasi-one dimensional compound family TiOX (X = CI, Br). With a simple model 
of spin- 1/2 chains inserted in a planar triangular geometry which couples them elastically, we are 
able to obtain the uniform-incommensurate and incommensurate-dimerized phase transitions seen 
in these compounds. Moreover, we follow the evolution of the wave-vector of the distortions with 
temperature inside the incommensurate phase. Finally, we predict gapless spin excitations for the 
intermediate phase of TiOX compounds along with incommensurate spin-spin correlations, offering 
to experimentalists a candidate to observe an exotic Luttinger liquid-like behavior. 

PACS numbers: 63.20.kk, 75.10.Pq, 75.40.Mg, 64.70.Rh 



I. INTRODUCTION 

The recent discovery of a new class of inorganic spin- 
Peierls (SP) materials, TiOX (X = CI, Br), has made 
reemerge the interest on this magneto-elastic transition. 
Unlike the previously known SP materials, these com- 
pounds present a novel intermediate temperature incom- 
mensurate phase between the usual uniform and dimer- 
ized ones. Other peculiar properties are the presence of 
a pseudo spin gap up to T* ^ libK— for TiOCl and the 
existence of a very large singlet-triplet energy gap in the 
dimerized phase with values Eg ^ A'iQK for TiOCl^ and 
2mK for TiOBr^. 

The magnetic Ti chains are arranged in a triangular 
lattice forming a bi-layer structure. Good Bonner-Fisher 
fittings to magnetic susceptibility measurements are ob- 
tained from which it is inferred that these compounds 
are well described by one-dimensional S — 1/2 Heisen- 
berg chains. Moreover, the exchange interaction con- 
stants which emerge from those fittings are J ~ QQQK^ 
and - 370iv:l for TiOCl and TiOBr respectively. These 
isostructural compounds present two phase transitions: a 
uniform-incommensurate one at Tc2 ~ and ~ ASK^ 
and an incommensurate-dimerized one at Td ~ QIK and 

21 K respectively. The order of the high temperature 
transition has been a subject of extensive controversy al- 
though now there appear to be enough evidence that it is 
a second order transition^ whereas that at Td is of first 
order—. 

We have recently presented a simple microscopic model 
of independent Heisenberg chains immersed in a two- 
dimensional triangular elastic lattice''. This model pre- 
dicts the appearance of the uniform to incommensurate 
phase transition through the softening of a phonon near 
the Brillouin zone boundary. This points to the rele- 
vant role of the lattice geometry as a driving mechanism 
for the incommensuration. On the other hand, a precise 



characterization of the magnetism inside the incommen- 
surate phase is still lacking due to the difficulty of ob- 
taining large enough single crystals to perform inelastic 
neutron scattering experiments. Recently, the dynamical 
structure factor of TiOBr was obtained by time-of-flight 
neutron scattering on powder samples^. These measure- 
ments suggest the closing of the magnetic gap when go- 
ing into the incommensurate phase from the dimerized 
phase. This seems to be in agreement with our theo- 
retical predictions which will be detailed in the present 
work. 

In the first part of this article we do a review of the 
mechanism leading to the incommensurate transition. 
The model of Ref. is extended to account for incom- 
mensuration in the direction perpendicular to the mag- 
netic chains as well. Within a model of XY spins coupled 
to adiabatic phonons, we obtain lattice configurations for 
different temperatures in the incommensurate phase. In 
the second part, we aim to characterize the magnetic ex- 
citations in this phase. We search for the existence or 
absence of spin gap and calculate spin-spin correlations 
functions with analytic and numerical techniques. 



II. LATTICE GEOMETRY AS A DRIVING 
FORCE FOR THE INCOMMENSURATION 

It is well known that a magnetic chain in which the 
spins are coupled to the lattice is unstable toward a 
dimerization of the ions. This is called a spin-Peierls 
transition. For a real compound to exhibit this kind of or- 
dering, quasi-one dimensionality is needed because large 
interchain magnetic interactions would make the system 
develop magnetic ordering before going into the singlet 
ground state. This is a possible reason why few SP ma- 
terials are known. Before the discovery of this transition 
in the TiOX family in 2003'^, the only known inorganic 
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SP compound was CuGeOs. This compound presents a 
magnetoelastic transition from the uniform phase at high 
temperatures to a dimerized phase as it is predicted for 
a single chain. So the question is the following: what 
makes TiOX behave differently so as to develop an in- 
commensurate phase in between? 

Let us consider the spring-like lattice model of Fig. [T] 
as representative of the Ti ions in the bilayer structure 
of TiOX. 



a) 




b) 




FIG. 2. Illustrative representation of the action of a dimeriza- 
tion on the interchain springs in the cases of a) square and b) 
triangular geometries. Dashed lines represent the interchain 
spring configurations after the distortion takes place. 



FIG. 1. Schematic representation of our simplified model. 
Only Ti atoms are included over the xy plane. K\, K2 and 
are the harmonic force constants acting when two atoms 
in the same chain, in neighboring chains and in next-nearest 
neighbor chains respectively move from their equilibrium po- 
sition, ai and a2 are the base vectors. The coordinate of an 
atom is given by Ri j = ia.\ 4- j&2 ■ 



normal coordinates. We want to study what happens 
with the elastic energy as we move away from a in-phase 
dimerization, i.e. from the component q = (0,7r). The 
square of the lowest dispersion curve'' for our model of 
Fig. [1] (see the phononic Hamiltonian Hp^ in Eq. [3]) 
reads 



We are considering that there are magnetic chains in 
the y direction which are only elastically coupled between 
them. represents the displacement of the i-th ion on 
the j-th chain from its equilibrium position and Xi, K2 
and K'i are elastic constants which act longitudinally. 
For simplicity, suppose that the ions move only in the 
direction of the magnetic chains. Fig. [2] clearly shows a 
marked difference between the situation where the ions 
in neighboring chains are aligned normally and that in 
which they are shifted by half lattice constant. 

If the ions dimerize in an in-phase pattern in the square 
lattice as shown in Fig. [2^, there is no cost of elastic in- 
terchain energy because the distance to the correspond- 
ing ion in the next chain remains unchanged. Therefore, 
the system will prefer this kind of distortion to lower its 
overall free energy. This is what happens in the case of 
CuGeOa. On the other hand, for the triangular geome- 
try (Fig. [^b) one sees that there is always an increase of 
interchain energy due to the alternating contraction and 
expansion of the interchain ionic distances. Thus, it is 
less favorable for the triangular system to dimerize as it 
has to pay more elastic energy than in the square case. 

Let us perform a more precise characterization of what 
would happen with the distortions in the triangular case. 
Now we are going to allow for movements also in the x 
direction (perpendicular to the magnetic chains) for more 
generality. The total mean energy of the elastic system 
can be calculated as Eph = | X]q ^^(q)!*?!*!)!^ where 
ri(q) is the phonon dispersion relation and Q(q) are the 
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where = 4Ka/M. If we develop this relation near 
q = (0, tt) we get 
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where S = n — qy. 

Therefore we see that the development has a linear 
term in S (or qy) in contrast to what happens for the 
square lattice in which the first non vanishing qy term is 
of second order. To first order, the square lattice does not 
alter its elastic energy in the vicinity of (5 = but there 
is an energy gain for the triangular lattice as long as the 
condition q^ < 2\/2 is fulfilled. This is not a very restric- 
tive condition as we are looking for solutions with qx ~ 0, 
i.e. a smooth interchain modulation. As qy is the com- 
ponent which couples to magnetism in our model, then 
if the magnetic free energy does not grow very rapidly 
with S we expect the total free energy to have a mini- 
mum shifted from 6 = 0. On the other hand, the first 
non vanishing q^ term appears at second order. Anyway, 
for a fixed ^ > 0, the system would gain elastic energy 
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is satisfied. In conclusion, this model would enable us 
to find an incommensuration in both g-directions as seen 
experimentally. 

For a quantitative characterization of the uniform- 
incommensurate transition, let us extend our model of 
Ref. m allowing ion displacements in the x direction as 
well. We introduce the following Hamiltonian, 



H — Hph + Hs + Hsph, 
2m 
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where i-^ j is the momentum of the i-th ion on the j-th 
chain, Ki^ K2 and are the elastic constants shown in 
Fig. [T]and a is the spin-phonon coupling constant. The 
spin-phonon part has been linearized, so the magnetism 
couples only to the y component of the ion movements 
to this order. It is useful to introduce a dimensionless 
spin-phonon constant through 



A = 



4a V 
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The phononic part of this Hamiltonian allowed us to cal- 
culate the dispersion relation ([T]). Following the steps 
of Ref. 1^ we calculate the phononic retarded Green's 
function. The phononic self-energy is treated by a Ran- 
dom phase approximation approach and we use the ex- 
pression obtained in Ref. HI for the dimer-dimer correla- 
tion function. There will be a structural transition when 
I?'''^*(q, oj) has a pole with w = 0. This gives us an equa- 
tion for the transition temperature as a function of q. 
The q which has the highest temperature will be the 
wave vector associated to the lattice distortion at the 
transition point. Fig. |3] shows a contour map of the 
frequencies of the softened phononic dispersion surface 
at the transition temperature, in which it is seen that 
the mode which drives the structural transition is incom- 
mensurate for both q-components, in agreement with the 
experimental measurements^^^'^. 

The free parameters i^i,i^2,^3,A of out model have 
been adjusted in order to fit different experimental data 
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FIG. 3. Contour plot of the frequency surface of the renor- 
malized phonons at the transition temperature in the area 
around q = (0, tt). The darkest zone encloses the wave vector 
which acquires the lattice modulation at the transition point. 



measured in TiOCl (cf. Fig. 3 of Ref. |6|): the phonon 
frequencies obtained by inelastic x-ray scattering for the 
(0, gy) direction at T = 300X^, the value of the bare fre- 
quency of the phonon identified as the one which drives 
the SP transition^, fio ~ 27meV, and {qx,qy) values sim- 
ilar to those found experimentally^iiSiM. Speci fically, t he 
values which gave us the best fit were: Vti = AKi /M ~ 
6meV, K2/K1 ~ 6.7, K3/K1 ~ 11 and A = 1.3. It is 
worth saying that this approach overestimates the tran- 
sition temperature by a factor ~ 2.5 as discussed in Ref. 
S We treated this problem in Ref. [l3 showing that an 
approach using non-adiabatic phonons solves the prob- 
lem, but we will not do the calculations for the present 
dispersion relation since it does not lead to qualitatively 
different results. 

Now let us turn to follow the evolution of the distor- 
tions for different temperatures inside the incommensu- 
rate phase. For this purpose we propose as a model an 
array of XY spin chains immersed on the elastic struc- 
ture of Fig. [1] As the XY spin-phonon chain contains 
the relevant physics of a SP system, we expect that the 
qualitative behavior will be captured by this simplified 
model, which incidentally can be solved exactly. The 
Hamiltonian is the one defined in ^ but now we ne- 
glect the z-component of the spin operators. Now we aim 
to minimize the total free energy (magnetic -|- elastic) 
self-consistently in the displacements at different tem- 
peratures. For this purpose we first perform a Jordan- 
Wigner transformation on the spin operators of this XY 
Hamiltonian, which leads to a tight-binding-like model 
of spinless fermions coupled to the lattice. This model 
can be solved exactly by diagonalizing Nc matrices of 
size Ns X TVs (where Ng , Nc are the number of sites and 
chains) each one, and whose elements depend on the ac- 
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tual displacement configuration. This allows us to obtain 
the full one-particle spectrum. Then, the free energy per 
site is obtained by the expression 

^ 1=1 m=l 

where Cph is the elastic energy density associated with 
Hph in Eq. [3] and A™ is the m-th eigenvalue of the ma- 
trix associated to the 1-th chain. By this minimization 
procedure we are able to obtain both phase transitions, 
i.e. we can account for the full intermediate incommen- 
surate phase. Fig. |3] shows two examples of distortions 
of one of the chains in the incommensurate phase, one 
near the dimerized phase and the other near the uniform 
one. Superimposed to each pattern, we show the wave 
which modulates the alternation (— together with 
a sinusoidal fit to it. 
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FIG. 4. Distortion patterns in the incommensurate phase, 
a) close to the uniform phase at T / J ^ 0.14 and b) near 
the dimerized phase at T/J ~ 0.12. Solid lines represent the 
waves that modulate the distortions given by (— l)'iti. Dashed 
lines are sinusoidal fits to the envelopes. 

We observe that in the case close to the dimerized 
phase, the modulation is not so well described by a sine 
function. The distortion pattern looks more like a soli- 
ton lattice which is confirmed by the secondary peaks 
appearing in the Fourier transform of the positions (see 
Fig. [51). 

When the temperature raises, the modulation evolves 
and close to the uniform phase it is a perfect sine function 
(Fig. [3^). This behavior could explain the experimen- 
tal temperature evolution of the x-ray peaks observed in 
Ref. ^9 for the incommensurate phase of TiOCl. There, 
a third harmonic of the main incommensurate peak was 
observed and interpreted as due to a soliton lattice struc- 
ture. Upon warming toward the uniform phase this third 
harmonic signal reduces and the pattern becomes sinu- 
soidal. On the other hand, this behavior is also similar to 
what it is observed in CuGeOa when a magnetic field is 
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FIG. 5. Fourier transform of the positions shown in Fig. [4)3. 
The appearance of secondary peak indicate the tendency to 
form a soliton lattice near the dimerized phase. 

applied to the compound in the low temperature phase^^. 
In that system, when the magnetic field is high enough 
to overcome the spin gap (first critical field), the system 
first develops a soliton lattice and with increasing field it 
tends to a sinusoidal pattern up to a second critical field, 
above which it becomes uniform. 

To finish this section, in Fig. [5] we show the evolution 
with temperature of the wave vector of the main Fourier 
peak. 
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FIG. 6. Evolution of S and Qx with temperature. Dashed lines 
are the frontiers which separate the different phases. 

We were able to obtain a similar dependency of qy with 
T as seen experimentally^iiS. On the other hand, the 
Qx modulation that we obtain increases slightly with de- 
creasing temperature, contrary to what it is found exper- 
imentally. This might be due to the fact that our model 
does not consider magnetic interchain interactions, so the 
system is not able to gain magnetic energy by a distortion 
in the x direction, in contrast to what occurs in the direc- 
tion of the magnetic chains. Fig. [7] shows the evolution 
of the amplitude of the distortion uq, in which we can ob- 
serve that above Tc2 we have uq = indicating that this 
is the order parameter for the uniform-incommensurate 
transition. 

The behavior of this curve around Tc2 resembles that 
of a second order phase transition (cf. with Ref. @). On 
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FIG. 7. Evolution of the distortion amplitude with tempera- 
ture. Note the finite jump of uo at Td and the second order- 
like behavior around Tc2- 
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FIG. 8. One-particle energy spectrum for the pseudofermion 
representation of the XY model with two different modula- 
tions. Both cases present continuous bands at the Fermi level, 
indicating gapless excitations. 



the other hand, at Td there is a finite jump of uq and 
q (Fig. [S]), signahng a first order transition. The order 
parameter for this transition is d. 

III. NATURE OF THE EXCITATIONS AND 
SPIN CORRELATIONS 

Let us now characterize the behavior of the magnetic 
excitations inside the incommensurate phase. We are go- 
ing to check for the existence or not of a spin gap and cal- 
culate spin-spin correlation functions. For this purpose 
we tackle the problem with two different approaches. 

In the previous section we have minimized the free en- 
ergy for a system of elastically coupled XY chains and 
obtained the corresponding distortions for different tem- 
peratures. Now we are going to use these configurations 
as stable patterns for one chain and study some prop- 
erties on it. First we will check the energy spectrum 
around the Fermi energy of the pseudo-fermion chains 
with hoppings modulated by these distortions. As dif- 
ferent Sz subspaces in the magnetic problem are related 
to a different filling of the one-particle energy levels in 
the spinless fermions, this will enable us to identify the 
presence or absence of a spin gap. In Fig. |S] we see the 
spectra for two different patterns. In both cases con- 
tinuous bands appear at the Fermi level, indicating the 
absence of a gap. We were able to check this behavior 
for several different patterns which always led to gapless 
spin excitations. 

Now we turn to the calculation of spin-spin correlations 
(SfSf^^) — (SfSf^j,). This is accomphshed through the 
calculation of r x r determinants whose elements are ob- 
tained by the formula— 

N/2 

G,,, = Yl (^"^r - er+^/'e7+^/') (9) 

m— 1 

where e™ is the i-th component of the eigenvector cor- 



responding to the m-th eigenvalue of the matrix intro- 
duced in the previous section. Fig. [HI shows real space 
correlations and Fig. [TU] shows the static structure factor 
along with the Fourier transform of the corresponding 
distortion pattern. 
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FIG. 9. Real space spin-spin correlation functions obtained 
with the XY model for the uniform chain and the distortions 
shown in Fig. ID 

On one hand we see that the spin correlations in the 
incommensurate phase are modulated. This is made clear 
by the appearance of new peaks in the structure factor. 
In every case there appear peaks at = tt — 2S and 
Qy — 6 (not shown in Fig. [TU]) besides the one at qy — tt 
which already belongs to the uniform antiferromagnetic 
lattice. With this information in hand we were able to 
fit the real space correlations by the following formula 

(5f5f^.) = ,^22^+,^ (-^y(^ + -^-os2Sr) ^^^^ 

where the parameters a^, which measure the relative 
intensity of each peak, take on different values for each 
distortion. On the other hand, the magnetic Fourier peak 
near tt is at twice the distance from tt of the lattice peak. 
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FIG. 10. Fourier transform of the correlations shown in Fig. 

[9] (solid lines) along with the peaks of the lattice modulation FIG. 11. Real space spin-spin correlation functions obtained 
(dashed lines). with DMRG for two different single-chain configurations (sig- 

naled with squares and circles) and the result for the uniform 
chain as a reference. 



It is worth mentioning that the peak at Qy = S has a much 
smaller intensity than those at = tt or (7^ = tt — 26. 

We now turn to the study of the Heisenberg model with 
modulated exchange on a single chain. Nowadays, the 
Density Matrix Renormalization Group (DMRG) numer- 
ical technique appears as the ineluctable numerical tool 
for treating static properties in one dimensional interact- 
ing systems at zero temperatur o^^'^^ . For our purpose 
we chose its finite size algorithm with periodic bound- 
ary conditions (PBC) on systems of sizes up to L = 200, 
instead of open boundary conditions to avoid edge ef- 
fects. Even though it is known that using PBC leads 
to a slower convergence of the method^^, keeping just 
300 states was enough to get quasi-exact results with 
a discarded weight of order 0(10^^^). In particular 
we are interested in the spin gap, here determined as 
A, = EoiS, = l)-Eo{S, = 0), being Eo{S,) the ground 
state energy on a chain in the subspace of a given total 
Sz ■ Different single-chain distortion patterns of the type 
Ui — uq cos(q- Ri) were tried and the absence of spin gap 
was found for every configuration, confirming a spin gap- 
less situation in the incommensurate phase. Spin-spin 
correlation functions have been also calculated and the 
results are shown in Figs. [TT] and[T2l 

As for the case of the XY model, the real space cor- 
relations show a power law behavior, indicative of the 
absence of spin gap. The decay is modulated by an en- 
velope whose Fourier component is q — t: — 26 as found 
previously. We propose, as a check for this prediction, 
the measurement of the static magnetic structure fac- 
tor by elastic neutron scattering. This should present 
an incommensurate peak at twice the distance from the 
commensurate peak of the low temperature phase with 
respect to that seen in the lattice structure factor. 

To conclude this section, let us comment on the results 
given by bosonization on the Heisenberg Hamiltonian and 
a Renormalization Group analysis. This procedure has 
been treated in Ref. [l^ for a problem of fermions with 
incommensuration. Here, we reproduce the bosonized 
version of the Hamiltonian using our parameters, which 




FIG. 12. Fourier transform of the magnetic correlations of 
Fig. [11] (solid lines) along with the peaks of the lattice mod- 
ulation (dashed lines). The symbols are in accordance with 
those of that figure. 



reads 



H = — dx 
2tt 



K 



2'Ka 



da; [sin(2(/) + (5a;) + sin(2(/) - (5a;)] (11) 



where (^ is a smooth bosonic field, H is the field conju- 
gated to 0, u is a velocity, K is the Luttinger parameter 
and a is a cutoff of the theory. 

The flow equation for 6 is given by ^ = (5 (cf. Eq. 



45 of Ref. l20f ). It is easy to observe that 6 grows with 
scale, which means that the sine functions in (|lip oscil- 
late rapidly so the integral will vanish. Therefore we end 
up with a free bosonic theory indicating that the exci- 
tations will be gapless as in a uniform Heisenberg chain. 
Note that the bosonized Hamiltonian of Eq. (|lip could 
correspond to a general XXZ model. The only change 
is the value of the Luttinger parameter K—. Therefore, 
the results presented in this Section seem to indicate that 
the behavior of the correlation function given in Eq. (jlOp 
has an universal form which is obtained by changing the 
exponent of the power law decay from a square root in 
the XX model, to a K-dependent exponent in a general 
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anisotropic XXZ model. 

IV. CONCLUSIONS 

In conclusion, we have studied the behavior of a sys- 
tem of spin-1/2 magnetic chains which are immersed in 
a triangular planar lattice. It is seen that this triangular 
geometry leads the system to a quite different behavior 
compared to a regular spin-Peierls mechanism, like that 
found in CuGeOs and other organic compounds. The 
elastic frustration produced by the triangular arrange- 
ment of the magnetic ions leads the system to produce an 
incommensurate phase between the uniform and dimer- 
ized ones. This is what it is observed experimentally 
in the compound family TiOX (X = CI, Br). We have 
seen that with a simple model of magnetically uncou- 
pled chains we were able to obtain both phase transi- 
tions. Moreover, by a self-consistent minimization of the 
free energy, we could describe the lattice configurations 
at different temperatures which allowed us to follow the 
y-componcnt of the main modulation wave-vector with 
temperature. On the other hand, spin-spin correlation 
functions were calculated exactly with XY chains and by 
DMRG calculations on the Heisenberg model on single 
chains with sinusoidal modulations. In both approaches 
we obtained power law behaviors at long distances, mod- 
ulated by an envelope which has a Fourier peak at 25, 
where at 5 is the peak of the imposed lattice modula- 
tion. Moreover, this power law decay of the correlations 



is indicative of the absence of spin gap in the excita- 
tion spectrum, which was confirmed by looking at the 
one-particle spectrum in the XY case, and the vanish- 
ing energy difference between Sz — Q and Sz — 1 sub- 
spaces in the case of the DMRG calculations. Further- 
more, an analysis of the bosonized version of the spin- 
phonon Hamiltonian together with the Renormalization 
Group flow equation of the incommensuration parameter 
5 led to the same conclusion. In recent time-of-flight neu- 
tron scattering experiments in TiOBr, it has been shown 
that the spectral weight increases at low energy when 
going from the dimerized to the incommensurate phas(^. 
This is in agreement with the theoretical prediction of the 
present paper that the gap closes in the incommensurate 
phase. A more detailed experimental determination of 
the magnetic spectrum and the theoretical calculation 
of the dynamical correlation function in the incommen- 
surate phase would be needed to further clarify the na- 
ture of the magnetic excitations in this uncommon phase 
which we have characterized as an incommensurate spin 
Luttinger liquid phase. 



ACKNOWLEDGMENTS 

The authors would like to acknowledge J. P. Clancy for 
helpful discussions. This work was partially supported 
by PIP 11220090100392 of CONICET, PICT 1647 and 
PICT R 1776 of ANPCyT. 



T. Imai and F. C. Chou, "Novel Spin-Gap Behav- 
ior in Layered S=l/2 Quantum Spin System TiOCl," 
'arXiv:cond-mat /0301425vl 

J. P. Clancy, B. D. Gaulin, C. P. Adams, G. E. Granroth, 
A. L Kolesnikov, T. E. Sherline, and F. C. Chou, "Singlet- 
triplet excitations in the unconventional spin-Peierls sys- 
tem TiOBr," arXiv:1009.4 185vll 
' A. Seidel, C. A. Marianetti, F. C. Chou, G. Ceder, and 



P. A. Lee, |Phys. Rev. B 6 7, 020405 (2003) 
R. Riickamp, J. Baier, M. Kriener, M. W. Haverkort, 
T. Lorenz, G. S. Uhrig, L. Jongen, A. MoUer, G. Meyer, 
and M. Griininger, Phys. Rev. Lett. 95, 097203 (2005) 
J. P. Clancy, B. D. Gaulin, and F. C. Chou, 



Phys. Rev. B 81, 024411 (2010) 



and 



A. 



Dobry, 



' Mastrogiuseppe 

|Phys. Rey^ B 79, 134430J 2 009)^ | 

Note that there are two acoustic dispersion branches within 
this model as we are working in two dimensions. We con- 
sider the lowest as the one which couples to magnetism 
more relevantly. 

M. C. Cross an d D. S. Fisher, 
|Phys. Rev. B 19, 402 (1979)] 

' E. T. Abel, K. Matan, F. C. Chou, E. D. Isaacs, 
D. E. Moncton, H. Sinn, A. Alatas, and Y. S. Lee, 
[Pi^^s. Rev. B 76, 214304 (2007) , 



A. Schonleber, S. van Smaalen, and L. c. v. Palatinus, 
Phys. Rev. B 73, 214410 (2006) 

van Smaalen, L. Palatinus, and A. Schonleber, 
iPhys. Re v. B 72, 020105 (2005) 

T. " Sasaki, M. Mizumaki, T. Nagai, T. Asaka, 
K. Kato, M. Takata, Y. Matsui, and J. Akimitsu, 
Physica B: Condensed Matter 378-38 0, 1066 (2006) 
proceedings of the International Conference on Strongly 
Correlated Electron Systems - SCES 2005. 
A. Krimmel, J. Strempfer, B. Bohnenbuck, B. Keimer, 
M. Hoinkis, M. Klemm, S. Horn, A. Loidl, 
M. Sing, R. Claessen, and M. v. Zimmermann, 
Phys. R ev. B 73, 172413 (2006) 

" " A. 



D 



Mastrogiuseppe 



and 



Dobry, 



jJournal of Ph ysics: Co nference S eries 200, 022032 (2010)^, 
I. Loa, The spin-Peierls system CuGeOi. A Raman- 
scattering study, Ph.D. thesis, Technische Univ., Berlin 
(1999). 
E. Lieb, 



T. Schultz, 



and D. 



Mattis, 



Annals of Phy sics 16, 407 (1961) 



U. SchoUwock. ^Rev. Mod. Phys. 77, 259 (2005) 
K. A. Hallberg, Advances in Physics 55, 4 77 (2006) 
S. R. White, Phys. Rev. B 48, 10345 (19937P~^ 
D. Sen and S. Lai, Phys. Rev. B 61, 9001 (2000) 
T. Giamarchi, Quantum Physics in One Dimension (Ox- 
ford University Press, 2004). 



